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_1.  _Introd':>;i.  i  on 

For  satisfactory  analytical  assessment  of  metal  forming  problems, 
it  is  necessary  to  evaluate  t)ie  varying  stress  distribution  throughout 
the  work-piece  since  1  arming  defects,  such  as  the  development  of  inter¬ 
nal  cracks,  depend  ori  .‘.tress  history,  and  residual  stresses  can  be  ii: ■>  r 
taut  in  deciding  the  utilization  of  a  formed  part.  Elastic  character i s~ 
tics  usually  play  an  essential  role  in  the  determination  of  stress,  oven 
in  combination  with  extensive  plastic,  flow  which  may  involve  strains  r. 
thousand  times;  elastic  strain  magnitudes.  Thus  analysis  of  metal 
forming  problems  for  .such  assessments  must  be  based  on  elaslic-plast j c 
theory.  The  same  is  true  for  other  stress  analysis  problems  when  plas¬ 
tic  flow  occurs,  unless  plastic  flow  is  occurring  simultaneously 
throughout  ilie  entire*  body  throughout  the  duration  of  the  process,  in 
which  case,  plast  ic-rigid  analysis  is  adequate. 

Because  plasticity  laws  are.  incremental  in  nature,  they  result  in 
relations  between  stress-rate  and  strain-rate,  or  equivalently  in  manor i 
cal  evaluations,  betas. on  stress  and  strain  increments.  For  the  rate- 
independent  laws  usually  adequate  at  temperatures  low  compared  with  th“ 
melting  temperature,  linear  relations  bet  •ecu  stress-rate  and  strain- 
rate  arise.  When  plastic  How  is  taking  place,  the  coefficients  are 
functions  of  the  current  stress  for  the  common  laws  and  below  tin*  yield 


stress  the.  elastic  laws  apply  in  incremental  form.  Because  of  the  struc¬ 
ture  of  t  hese  laws,  elastic-  pi  ast.ic  problems  are  commonly  solved  in  terms 
of  equations  for  stress-rates  and  strain-rates,  containing  stresses  as 
coefficient s.  Consider  the  solution  to  have  been  computed  up  to  the 

time  t  .  After  a  time  step  forward,  At  ,  the  solution  for  rates  gives 

» 

the  stress  at  time  t  +  At  in  the  form  o(t)  +  oAt  ,  and  similarly  for 

other  variables.  0  i;i  *  he  appropriate  stress-rate.  Then  a  new  tine 

step  can  le  to  hen  end  the  process  repeated. 

Extensive  studies  oi  Lho  application  of  such  laws  to  stability  and 

uniqueness  of  solutions,  have  been  made  by  Hill  (see,  for  example  [1),  [7.] 
*> 

[3]  where  he  shows  tint  care  in  the  selection  of  stress  definitions 
and  stress-rate  and  strain-rate  expressions  is  important  for  a  satis¬ 
factory  development  of  the  theory.  Rice  [4]  has  pointed  out  that  such 
questions  are  also  important  in  developing  a  satisfactory  theoretical 
basis  ior  elastic-plastic  stress  analysis,  particularly  in  the  common 
circumstance  that  the  tangent  modulus  in  plastic  flow  is  of  the  order 
of  tho  stress.  Converted  and  rotation  terms  then  become  important  in 
the  stress  rat e  expresr’on,  and  analogously  stress  variables  should  be 
selected  so  that  the  Influence  of  rate  of  deformation  of  the  boundaries 
of  the  body  does  not  complicate  the  variational  principle  commonly  used 
to  replace  me  oquil ibi ium  equations  for  evaluation  of  solutions.  This 
requirement,  can  he  achieved  by  using  the  unsymnctric  nominal  stress 

k 

numbers  in  square  brachets  refer  lo  the  bibliography. 
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(Piola-Ki rchof f  I)  in  which  the  r.trcss  is  defined  as  force  per  unit 


undeformcd  area.  The  variational  principle  then  involves  an  integral 
over  the  uridcformed  body  which  is  fixed. 

Plastic  flow  is  essentially  a  fluid  type  phenomenon  which  can  be 
most  conveniently  expressed  in  terms  of  the  current  configuration  of  the 
material.  Thus  a  reference  configuration  which  remains  invariant  through¬ 
out  the  motion  is  not  appropriate  and  so  the  configuration  at  time  t 
is  adopted  as  the  reference  ?  for  evaluation  of  the  deformation  from 

t  to  t  +  At  ,  where  At  is  efficiently  small  for  adequacy  of  first 
order  theory. 

The  framework  described  above  provider,  n  satisfactory  foundation 
for  a  finite-element  elastic-plastic  code  as  discussed  by  McMeeking  and 
Rice  [5].  In  effect,  by  choosing  the  current  configuration  as  the  refer¬ 
ence  state  the  Cauchy  stress  (or  true  stress  in  Cartesian  coordinates) 
the  unsymmotric  nominal  stress  (Lagrange  or  Pioln-lCirchhoff  I)  and  the 
symmetric  nominal  stress  (Kirchhoff  or  Pi.ola-Kirc.bof  £  II)  all  have  iden¬ 
tical  values  at  the  current  time  which  tjimplifi.es  utilization  of  the 
appropriate  stress  for  the  appropriate  component  of  Lhe  calculation. 
Although  the  stress  components  themselves  are  identical,  rates  ot  change 
of  the  different  stresses  are  not  the  same. 

2.  Dove! opment  of  the  Then ry 

Following  Hill  [?.]  and  using  for  the  most  part  his  notation,  we.  con¬ 


sider  the  unsymmetr ic  nominal  stress  (variously  referred  to  as  Lagrange 


or  Piola-KLrchof f  I)  s  cU-fined  so  that  the  jth  component  of  the 
force  transmitted  across  a  deformed  element,  which  in  the  initial  or 


reference  state  had  area  dS  and  unit  normal  ,  la 


dS  v  S. .  »  dr.  (1 

i  iJ  .1 

Hill  considers  (p.  214  of  (?.])  rate  or  flow  type  constitutive  la\;s  of 
the  form 


ij  ~  3(Dv  /9X.) 

where  X.  are  rectangular  Cartesian  coordinates  in  the  initial  or  refer 

once  configuration*  F,  is  a  homogeneous  function  of  degree  tv;o  in  the 

velocity  gradients,  9v  /9X.  >  and  where  s..  is  the  partial  time  deri- 

j  i  ij 

vativc.  at  fixed  X  ,  i.e.  a  material  derivative  at  a  particle.  The 
velocity  v. (X,t)  gives  the  distribution  at  time  t  expressed  in  the 
initial  coordinates  of  the  corresponding  material  points  (note  that  a 
tilde  under  a  symbol,  denotes  a  vector  or  tensor  in  absolute  notation). 

O 

Boundary  value  problems  are  considered  in  which  for  a  volume  V 
in  the  reference  state,  at  time  t  stress  rates  (X,t)  and  veloci- 

ties  v^(X,t)  are  sought  for  prescribed  nominal  traction  rates 

O 

over  the  part  of  the  vfnee  S„  ,  velocity  v.  over  the  remainder  of 

1  J 

«  . 

the  surface  and  body  force  rates  g^  per  unit  initial  volume. 

Then  the  variational  principle 


6 1  /  E  dV  —  /  F.v.dS  -  /  g.v.dVj  =  0 
o  o  J  j  0  j  j 


in  the  class  of  continuous  differentia bio  velocity  fields  satisfying 


the  velocity  boundary  condition  on  ,  characterizes  the  solution, 
for  it  yields;  the  equilibrium  equations; 


129. 


—±i 

ax. 


+  e.  ~  0 


(4) 


and  the  boundary  r  motion  rate  condition 


Vij 


for  nominal  stress,  s  ,  and  the  reference  geometry. 

In  writing  the  elastic-plastic  constitutive  relation,  we  wish  to 
associate  the  veluaity  gradient  in  (2)  with  the  .rate  of  dcfoi .nation  or 
velocity  strain: 


D. . 

ij 


3v.  3v, 

(v.-V 
«?:.  Dk/ 
J  i 


where  x^  are  Cartesian  coordinate:; 
cles  in  the.  di  formed  body 


expressing 


the  position  of 


(6) 

par  1 1- 


Xi  =  Xi(X;i:)  (?) 

Thus  to  permit  simultaneous  use  of  (3)  and  the  plasticity  laws  expressed 
in  the  usual  form  of  rate  of  deformation  of  the  current,  configuration. 
Hill  takes  the  current  configuration  to  be  the  reference  state,  and 
hence 


X.  -  x.(X,t) 
i  x  -  . 


(8) 


for  a  particular  t  .  The  theory  is  expressed  in  this  fern  for  evalua- 
tion  of  and  v  and  hence  the  solution  and  configuration  at 

t  +  At  ,  which  provider,  a  new  reference  state  for  evaluation  of  the 


next  tine  step  At  . 


Note  tin:  at  !  he  jnst.m: 


t 


when  the  current  and  reft  mure  con¬ 


figurations  arc  i.,i>  n.L ical ,  the  nominal  stress  components 


i.i 


are  equal 


to  the  Cauchy  or  u  ue  stress  components  >  so  that  at  this  instant 


The.  dev  ice  of  selecting  the  coni  igur.it ion  at  time  t  to  be  the 

reference  .  rate  e.va !  u.tr  !;  a  el  the  solution  at  time  t  +  At.  thus 

perm 1 1  simul  tain  ■■■.<,  use  of  the  convenient  variational  principle  (3)  in 

terms  of  nominal  .-.t  ress  s .  .  and  a  fixed  geometry  and  the  familiar 

plasticity  Jaws  e  pressed  in  terms  of  the  Cauchy, -or  true  stress,  . 

By  working  Lrh  curvilinear  convi  cted  coordinates,  ,  having  an 

arbitrary  configuration  in  the  refer.r.t,  or  initial  state,  Hill.  ([2], 

■  •% 

p.  21'j  if.)  shows  that  tli  rate  potential  (2)  follows  from  associated 
rate  potentials  for  other  .stresses  end  stress  rate  expressions,  some  of 
which  arc  more  convenient  for  representing  elastic-plastic  laws. 

Consider  cm  vi  linen  i  eo.ird  in.  t  <•:.  in  the  initial  or  reference 

C 

state  with  l>a-:u  tor*'  g  .  After  t!  -f  ormat  ion  a:,  shown  in  l'ig.  1 
these  become  C  ,  having  the  sane  values  for  the  same  material,  par¬ 
ticles  as  .  Then  are  the  convected  coordinates  and  the  cor- 

O 

refiponiing  h.n.e  vectors  are  g  ,  v/hlcii  are  g  deformed  by  the 

~u'  -u 

motion.  Then  the  Cauchy  stress  tcnsoi  a  has  contrravariant  components 

q  +  a  * 

o  in  converted  coordinates  such  that  the  force  dK  transmitted 

across  an  element  of  the  dr  formed  hotly  ol  area  tlS  and  unit  normal 

vector  v  ,  Is  given  by: 

a 


di'  •=  cQ  l'  (v  ,dS;g 

Cl  -i, 


(9) 


i  n . 


where  J  is  the  Jacobi  u  oi  the  transformation  from  the  reference  to 
the  deferred  state,  i.r. 


r  •  p  >:  p  3  j(o  *  v  x  o  ) 

"u  Ch  "y  is  -.'3  _y 


Cl/.) 


Fig.  2  shows  the  particular  situation  when  the  4°  coordinates  in 
the:  reference  frame  are  Cartesian,  X°  .  The  Cartesian  coordinates 
representing  points  in  the  deferred  body  .according  to  the  point  trans¬ 
formation,  e.-ju.  (7) ,  arc.  x  .  The  usual  definit  *on  of  the  Kirchhoff 
or  Piola-Ki rchho f f  II  stress  (s  o  Fung  [6]  p.  439)  is  given  in  terms  of 
this  point  transformation.  From  cqu.  (7),  define  the  deformation 
gradient 

i 


F-  12- 
3X° 


(15) 


Then  the  Piol a-K irchhof f  II  stress  r  is  given  in  terms  of  the  Cauchy 

c 

stress  in  Cartesian  coordinates  x  ,  a  by 


—  1  1  1 

T  “  d  F~J  a  F  ([7]  p.  125) 


(16) 


q  3 

where  J  ~  del:  (F)  .  "a is  is  in  accord  with  (13)  where  cr  are  the 
contravariant  components  of  the  Cauchy  stress  with  respect,  to  the  con- 
vcct.ed  coordinates  X  ”  for  if  a1"1  are  the  Cartesian  coordinates  of 


a  with  respect  to  x  ,  then  the  tensor  change  of  variabies  law  gives 

(17) 


a'fl'  ?,yu'  dxB  Ucij 


Sx1  3x1 


in  terms  of  the  coordinate  transformation  from  x  ■+  X'  in  the  deformed 
geometry.  Rccau.se  of  the  property  of  convectcd  coordinates  that  X°  =  X° 


.1  33. 


the  coordinate  transformation 


for  the  same  p-ut.’rlo,  (7)  also  expresses 
and  (13)  anu  (17)  are  seen  to  he  equivalent 
has  been  pointed  out  by  Nonet- Hamer  [8], 
tensor  densities  and  not  absolute  tensors, 
are  not  pure  tensor  relations. 


to  (16).  This  connection. 

v,  ,  <*6  ,  <:f> 

H o t: e  that  n  and  t  n re 

so  that  equations  such  as  (7) 


For  rate  inn. -pendent  constitutive  relations  the  rate-potential  lun  - 
tion  E  in  (2)  in  a  homogeneous  function  of  degree  two.  For  the  choice 
that  the  ro.feronn  >  state  is  the  current  state,  (2)  follows  from  the  exis¬ 
tence  of  an  associated  rate  potential  function  F(e  )  ,  of  the  strain  - 

op 

rate  component:; 


“op 


(v  +  v  ) 
a,B  8, a 


(th) 


where  v  is  the  velocity  vector  in  the  deforming  body  and  the  comma 

•  *  ct  fi 

denotes  covariant  differentiation.  This  function  generates  x  through 


•«R  _3F 

9  e 


c*S 


0  9) 


where  the  superposed  dot  indicates  time  derivative  of  the  convected  com¬ 
ponents  or  convected  derivative.  This  is  equivalent  to  the  partial  time 
derivative  at  f^xed  X  ,  or  a  material  type  derivative.  The  structure 
of  (19)  indicates  that  since  t  is  a  tensor  density,  F  is  a  scalar 
density  and  not  an  absolute  scalar  invariant. 

A  derivation  of  a  relation  needed  in  the  following  analyses  to  estab¬ 
lish  (2)  is  given  in  the  Appendix.  It  is  the  relationship  between 


•  ci  * t  *  Cl  p> 

s  *'  and  i  ,  which,  f  or  the  particular  choice  of  ref  or. -ace  state 


ment : oned ,  takes 


the  form 


decree,  rate-pel  n.!  -O 


s 


aft 


Loa  of  u,  for  tin-  nominal  stress  rata 

ft,  a 


A  similar  .  rgumenr.  f  s>y  the  Cauchy  stress  oC'^  yields  a  trial  pc 

tiel  function,  wi.;  ch  ho;. ever'  fa j  1c  to  yield  the  correct  expression  i'c 

•  ap 

o  from  the  rm  e-potenrl'.l  equation,  ller.ee  a  rate  potent  ial  fur.cti 
doer;  not  exist:  in  this  c a 

The  lav:!.,  of  plasti city  are  normally  obtained  by  measuring  "tru.; 

stress"  and  increments  of  attain  defined  in  terms  of  Cartesian  coord! 

nates  in  the  current  configuration  without  rotations  occurring.  Sine 

the  theory  must  apply  in  the  presence  of  rotations,  their  influence 

must  not  affect  the  stress  rate  term  in  the  constitutive  re.3r.ticn,  h. 

a  spin-invariant  stress  rate  is  needed,  such  as  the  Jaunanr.  rate.  We 

will  work  in  terms  of  the  configuration  at  time  t  ,  which  is  the  1 ef 

ence  configuration,  and  utilise  Cartesian  coordinates  x  . 

In  formulating  the  finite  element  theory  for  numerical  analysis 

elastic-plastic  problem*,  we  vish  to  use  the  variational  rotation  (3) 

in  terras  of  s..  arid  a  constitutive  relation  associated  with  (19)  in 
xj 

terms  of  t  since  it  can  be  conveniently  associated  with  measurerae 

of  plasticity  laws.  We  have  seen  that  a  law  in  the  form  (19)  implies 

the  validity  of  (2)  and  hence  the  variational  principle.  As  pointed 

out  in  (5]  this  structure  in  terms  of  t . .  leads  to  symmetric  stiff- 

rj 

ness  matrices  in  the  finite  element  formulation  which  simplifies  the 
numerical  procedures. 

I.’ow  the  rclaLioiw.hip  between  tin;  nominal  stress  s  a  id  the 

c 

Cartesian  true  stress  d  is 


lib. 


F  s  -  J  o 


see,  for  example  [7],  p.  3. '*‘5  where  the.  nominal  stress  is  defined  a;>  tin. 
transpose  of  t;  (or  :U  c.  e.  bo  deduced  from  (9),  (10)  and  (17)).  Taking 
the  material  do  rive  five  o.f  (24) ,  and  noting  that  F  and  J  are  unity 
for  coincident  re.feivwe  on!  current  configurations,-  one  obtains 
?v 

•  .  I  c.  k  ,  c  ,  „ , 

s.  .  +  s  .  »  o.  .  % —  -!•  o.  .  t?.3) 

13  kj  °;'j.  -J  "::k  3J 

The  difference,  between  the  Jaumann  derivative  of  d  and  its  material 
derivative  is  the  contribution  of  tin:  rate  of  rotation  of  the  axes  which 
rotate  with  the  body  according  to  the  nnti-syrmatric:  tensor  angular 
vel ocity  expression 


where  A  denotes  the  ant i-symmo trie  part.  Using  3)/i)t  to  denote  the 
Jaumann  derivative,  this  determines  (see  Prage.r  [9j  p.  lbd) 


S)o.  .  ♦ 

bl  c 

-J-  _o 


ik  ux, 


'kj  «)x. 


Combination  of  (23)  ,  (7.7)  and  (3) 


:»o .  .  dv,  ov, 

s,  ,  (  ---  --  +  S. ,  ---  )-  (a , ,  D  ,  +  a,  ,  D._  )  +  a.,  --J- 

ij  ‘J!t  xj  <»x,_  jk  jk  kj  ik  ik  ux, 


Equation  (13)  defines  t  j  ,  the  Kirchho.  f  stress,  in  terms  of  con- 
vected  coordinates,  l’.eing  a  tensor  density,  definition  for  other 
coordinator  is  obtained  by  use  of  (17)  suitably  modified  to  incorporate 


the  density  term  J  .  lie  have  seen  that  experiments  from  which  the  laws 


137. 


of  plasticity  were  i’w.o.:  involved  "true  sc  rear."  associated  wit  is 
Cartesian  coordinate-  which,  tor  plastic  analyses,  indicates  the  appro¬ 
priateness  of  a  Jan. -.am  tine-df  rlv.it  i  ve  associated  with  rotation  of 

rectangular  axes.  Thus  we.  need  to  utilize  r..  -  J  l'  to  incorporate 

xj  ij 

the  Jaurnann  derivative,  and  the  f  irst  ten:  in  parenther.is  in  (2d)  can  bo 
written  ‘J)T^./'i)t  ,  since  the  derivative  of  the  scalar  density  J  is 
unaffected  by  rotation  of  axes,  a, id  at  tine  t  the  instantaneous  value 
of  J  =  1  .  Hill  lias  stated  !  "J  p.  222)  that  the  rate  potent  fai  ( 1  2) 
implies  a  rate  potential  lor  the  Jaurnann  derivative  of  i  ,  which  pro¬ 
vides  computational  advantages  .associated  with  the  use  of  the  Juursnrt 
derivative  of  this  stress  variable. 

Thus  (28)  takes  the  form: 


:>Vi 

s  .  -  — — J  -  (c*  D. .  +  a,  .  D . ,  )  +  u  . ,  s  ** 
ij  !t)t  ik  jk  kj  x.<  ik  dx^ 

How  combining  (20)  and  (29) 


5)r .  . 
J_L1 
5)t 


Tij  +  °ik  Djk  +  0kjDik 


(?9) 


(30) 


In  view  of  (19)  written  for  initial  Cartesian  coordinates,  with 

F(D  )  homogeneous  of  second  order  in  the  strain  rates,  multiplication 

of  (30)  by  I).,  and  using  Kuler's  theorem  gives 

rf..  ,J 

S”1  D.  -  2F  +  o .  .0  D. .  +  o  D  D.  .  -  2F  +  2o  D.  .D  (31) 

ij  ij  jk  JJ  kj  ik  ij  ik  ij  jk  v  ' 

which  analogously  to  (21)  for  sa®  provider;  a  trial  rate-potential 
function 

G  -  F  +  c..  D. . D  (32) 

xk  ij  jk  v  f 
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assuming  Euler's  theorem.  In  a  manner  similar  to  that  presented  for 
s  '  ,  a  third  raie-potentla.1.  function  is  thus  established: 


!Li .  » 

•J)t 


(33) 


ij 


The  classical  elastic-plastic  isotropically  work  hardening  lav 
(Prendtl-Eeuss) ,  commonly  giving  strain-rate  as  a  linear  function  of 
stress-rate,  cm  be  inverted  to  give  stress-rate  and  takes  the  form 


<[b]  p.  606) 


_ 


sr  ■  #»•  iVn +  7. 


6  _  3ai;j  °k3  (  l  i'v  ) 


•2v  a  kl  2-2  ^  + 


1 


kl 


(34) 


where  E  and  v  are  Young's  modulus  and  Poisson's  ratio,  a'  denotes 
stress  devi.vcor,  a  is  the  current  tensile  yield  stress  and  h  the 
current  gradient  of  the  true-stress  logarithmic  plastic  strain  curve  in 
a  tension  test.  The  Jaemnnn  derivative  is  used  for  stress-rate  as  men¬ 
tioned  above  in  order  to  eliminate  rotation  effects,  and  the  last  term 
in  the  brackets  is  dropped  when  the  increment  of  deformation  is  elastic. 

c 

We  have  shown  that  there  is  a  rate-potential  function  G  for  t 


hur  not  one  for  o, 


i-J 


ij  ’ 

This  means  [5]  that  a  non-symmetric  finite- 


element  stiffness  matrix  would  be  deduced  using  (34).  Replacing 
c  c 

by  t  in  the  stress-rate  term  would  yield  a  symmetric  stiffness 
matrix  which  would  simpliiy  the  numerical  analysis.  Moreover  such  a 
change  Is  appropriate  in  terms  of  its  representation  of  the  physical 
laws.  The  J  term  in 


c  T  c 
T  .  .  =  J  0  .  . 
i.l  ij 


(35) 
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arises  in  wore  accur.o  v  i  ,  iresontai  ion  of  the  J avia  of  elasticity  than 
Hooke’s  law  used  in  c i.vu.Aenl  elasticity.  It  is  associated  with  gee— 
metrical  non-linearity  which  cxprenr.es  the  non-linear  influence  of 
finite  strain.  In  effect  it  e  >t v:;.xs  the  fact  that  energy  density  pec 
unit  initial  aolu.no.  yields  a  sv-.ph-r  energy  conservation  statement ,  lot 
per  unit  current  volume  ir.u>li?.s  change  of  energy  density  a  imply  because 
the  amount  of  material  containing  it;  changes.  For  example,  the  term 
(p/p  )  appearing  in  equation  (''6)  of  f  10)  is  equivalent;  to  replacing 
0  by  .To  ,  and  it  w.is  jn-.-.u J.  »•  jot  in  that  pnpe.j  (p.  935)  that  such  a 
term  provides  a  good  approxjji.tion  to  non-linear  elasticity  of  met; 1 s 
with  input  of  only  the  two  classical  elastic  constants.  A  similar 
modification  of  the  classical  laws  of  plasticity  was  suggested  in  [11] 
where  it  was  found  compelling  to  express  the  yield  stress  in  terms  of 
Jo  (eque.  (33)  and  (3-))  of  (11]).  Again  this  was  based  on  t.he  require¬ 
ments  of  geometrical  non-linearity,  which  of  course  are  independent  of 
specific  material  characteristics. 

Introduction  of  (2)  into  the  variational  principle  (3)  expresses 
it.  in  the  form 

i)V  ,o  #  o  o 

Is,.  6(v— L)  dV  -  &(f  F.v.dS  +  I  g.v.dV)  ■=  0  (36) 

o  rj  o  j  3  «  3  J 

V  SF  V 

which,  si  the  rcioj  once  state  is.  the  current  state  (X  ~  >:)  can  bo 

written 


Ov. 

I  s.  6  (~  *1')  d  V  -  <S(/  F.v.dS  +  /g.v.dV)  »  0 
V  3J  Dxi  s  3  J  v’3  3 

Substitution  of  (29)  and  applying  algebraic  manipulation  based  on 


(37) 


1  AO. 


symmetries  then  yields  tin?  vari.it  ioi.nl  principle  in  the  form  (a*  a  c  u. 


C>)  of  [3  1) 


/ 

V 


>'(D.  .) 


hP,  *(. 

Kl 


2D  ,  I, 
ik  kj 


k,  1 


v  . )  ]  dV 

k,J 


where  the  sub  :c  rl.pt 
(33),  the  principle 


-6 ( /  y.v  ds  +  /g.v . (IV) 

c  -1  J  o  J  J 


(38) 


,i  denoter.  the  operation  3/3x^  .  Now  utilising 
takes  the  form  (see  c.c;u.  (13)  of  {3]). 


<5  (  /  f-(h;  tlv  -  /a,  .  (2D  h  , 


.  v,  ,)d</ 


ij  ik  kj  k,i  k,j 


-  /  F.v.dS  -  / R.v.dV)  =  0 
S  J  J  V  3  J 


(39) 


My  (33),  am!  i; sing  Euler’s  theorem 


3?C 


c 

SIT.  , 

— =J-  =  ^  — ~_r - D  -  i’  D 

5)t;  an, .  3u..ah,1  ki  ijki  n 

ij  ij  kl  J 


(AO) 


since  3C/ai)_  is  homogeneous  of  first  degree  in  D  ,  hence  .9^  ,  ^  is 
synunetric  in  ij  kl  as  well  as  i  ■*->  j  and  k  «-*■  1  .  Since  G 


is  homogeneous  of  second  degree 

5)T., 

C(W  -  1;  Dy 


(41) 


The  symmetry  implicit  in  the  variational  principles  (38)  and  (39)  with 
(40)  and  (4i)  Imply  symmetric,  stiffness  matrices  which  carry  tin  c  e;  h 
to  the  finite-element  formulation  [5]. 

3.  .  Dit.cur.sicn 

As  discussed  by  Rice  [4]  and  McMeeking  and  Rice  [5],  the  develop¬ 
ment  reviewed  in  these  notes  brings  out  the  importance  of  convection 
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1 


effect!:  and  appropriate  ot  t .  .o-rate  definitions  in  i  oiruili.e  ;.in,  cln.lJc- 
plaattc  theory.  Tim:,  the  precise  analysis  of  cni.bLiiuont  the' : ^  j.  .ra  W 
applied  to  ohm  I  tel  li.Mc  r. suits.  Co:r..f»ji  small  strain  nr;:. a.  .prion  •, 
arc  not  vo  1  id  u  a  for  lari  octal  theory  based  on  small  strain  J.nci  u- 

rent:..  Tl-.ls  :  f  i  :.*«tLcn  is  1  licit  in  comparison  of  tl.c  first  and 
second  t  err.::  in  the  first  volume  integral  in  the  variational  principle 
(38) .  in  plan!  :c  r Jott  the  coefficient  of  6D  in  the  first  term  is 
hi)  ,  where  ;«  is  the  fp.iuient  of  the  tensile  stress-plastic  strain 
curve,  •  -I.  In  second  t<  v.i  is  *  o'lijf)  „  The  second  tern  tea  the 

diffeit-  ce  hefv.-en  the  Jaui;i.  mn  derivative  and  other  simpler  tfite.  der- 

vatives  can  o.‘y  be  neglected  if  h  »  o  .  The  relative  error  in 

neglecting  noth  terms  is  effectively  independent  of  the  magnitude,  of 

the  .strain  i ncrei lent  adopted,  su  that  small  steps  do  not  permit  uir, 3 i 

fleet  ion  In  thl  i  regard.  For  many  metals  lr  -  o  . 

It  is  interesting  to  note  that  the  complied  ions  which  arise  in 

elastic-plastic  analysis  result  from  the  elastic  component  of  strain, 

and  not  the  plastic.  For  stress  and  strain  deviators,  clast Ic.-icleatly 

2 

plastic  deforr  it  ion  with  a  Micas  yield  condition,  J  --  o'  a'  }7  «  1 

1  ij  ij 

satisfies 


’  “Ij/2C  +  Sj 


where  X  is  a  parameter.  Multiplication  of  (42)  by  o'  gives 

•  *  J 


(«; 


x  -  oijnj'j/2J2 


(43) 


For  isotropic  work  hardening  with  o  Mises  yield  condition,  (42)  takes 
the  form 


1/.?. 


.  .  ,.8*; 


"  /:>C.  ;•  1  {j 

•hi 


(44) 


The  first  It..-;  rn  t>.  right  hand  .sides  of  (42)  and  (44)  arc  Che  c!i.  Stic 
strain  rote  component:  .  *1  iu.s  rigid-ideally  plastic  theory  ((42)  w.t.h 

the  u  torn  .i.Jctrd)  ,-,'vvs  a  relat  ion  between  stress  .arid  velocity  gra¬ 
dient  with  no  pi. i.  o. 'ons  due  to  a  stress  rate  term,  which  grc.  i  ly 
simplifies  the  r  aaly  !  apart,  iron  the  difficulty  of  determinin',  the 
rip, id  iif.idns.  S  -  ai '•  Vv  for  work  hardening  rigid-p?  s'.ic  aua'iy 


((44)  with  tin:  o  te.ru  dela  ted),  only  the  rate  of  cl  ••  nye  of  a  ; 

invariant  occurs,  which  is  simpler  t.o  include  than  a  tensor  rat"..  .'or 

elastic-plastic  theory  j.v.  is  the  elastic  term  which  introduces  sir 

rate  and  the  consequent  complications. 

Many  technological  l.y  important  metal-  forming  problems  are  sf.  .  ciy 

state  proccssa  s  in  which  So,  ./l*t|  -  0  .  in  planning  to  use  an  as.  ar¬ 

id  1  x 

lysis  of  the  type  considered  here  to  evaluate  such  situations,  it  is 

fortunate  that  the  stress-rate  terms  appearing  in  the  variatior.nl 

•  c 

principle,  ,  or  3)i \^/£>t  ,  do  not  approach  zero  in  the  steady  case, 

hence  singular  computational  conditions  need  not  bo  anticipated.  This 
is  not  the  case,  in  some  simpler  and  inadequate  appro  aches  to  t.lrls  prob¬ 
lem  in  which  sufficient,  care  was  not  devoted  to  the.  appropriate  choice 
of  stress-rate  definition. 


On  the  basis  of  the  theory  described,  a  finite-element  program 
has  been  written  to  evaluate  stress  and  deformation  histories  in  an 
extrusion  problem.  The  case  of  a  billet  being  pushed  through  a  die 
until  a  steady-state  configuration  was  reached  has  been  completed. 
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The  Stic:;?;  field  exhibits  i  (Mtuces  which  arc  cap:;  j  tent  with  l  he  i  r.own 
development  of  extrusion  «t • : *  <  ets,  such  as  the  appearance  of  so r l  >■  c 
cracks . 

4.  Ach.iowl  cdgo.in.ii! 

In  study  in}’  (a-  n.nteriaj  preparatory  to  writing,  this  rev iVw,  i 
benefited  from  helpful  c.onvei.\  .ions  with  Professors  J.  W.  Hutch i r  ms , 

R.  .1.  Kaul,  R.  ],.  J-J.il  l.ett ,  J.  ;i  vers  and  S.  Neinat-Uasser  of  Harvard, 
State  of  Row  York,  St.:nlr*-d,  :  •  riford  and  Rorthv.-cst  era  Universal  > 
respectively,  and  with  Dr.  it.  ...  McKeekinj,;  of  blown  University 

This  work  was  sponsored  ia  part  by  the  Materia!:;  Science  PfvtsJ 
U.  S.  Office  of  Naval  Research  on  a  contract  with  Stanford  University, 
and  was  carried  out  in  part  while  the  author  held  a  dohn  Simon  Cu/.gi  nhei 
Memorial  Fellowship. 
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6.  Appe  i  u!  i  x 

For  simplicity,  relations  are  developed  for  Cartesian  arcs  in 
the  reference  state  which  is  instantaneously  coincident:  v;lth  the 
current  state.  The  i’.eaeral  theory  for  eonvoctcd  coordinates  carries 
through  in  a  similar  manner  but  can  be  technically  much  more  involved. 
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